Dynamics of Skyrmion Crystal in Metallic Thin Films 
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We study the collective dynamics ol the Skyrmion crystal (SkX) in thin films of ferromagnetic 
metals resulting from the nontrivial Skyrmion topology. It is shown that the current-driven motion 
of the crystal reduces the topological Hall effect and the Skyrmion trajectories bend away from the 
direction of the electric current (the Skyrmion Hall effect). We find a new dissipation mechanism 
in non-collinear spin textures that can lead to a much faster spin relaxation than Gilbert damping, 
calculate the dispersion of phonons in the SkX, and discuss effects of impurity pinning of Skyrmions. 
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Introduction: Skyrmion is a topologically nontrivial 
soliton solution of the nonlinear sigma model. It was 
noted early on that Skyrmions in three spatial dimen- 
sions have physical properties of baryons and the peri- 
odic Skyrmion crystal (SkX) configurations were used to 
model nuclear matter [l|, Skyrmions in two spatial 
dimensions play an important role in condensed matter 
systems, such as quantum Hall ferromagnetsd, |3|- It 
was suggested that SkX configurations can be stabilized 
by the Dzyaloshinskii-Moriya (DM) interaction in ferro- 
magnets without inversion symmetry |5| . Such a state was 
recently observed in a neutron scattering experiment in 
the A-phase of the ferromagnetic metal MnSi|6]. 

Recent Monte Carlo simulations indicated much 
greater stability of the SkX when a bulk ferromagnet is 
replaced by a thin film [7]. This result was corroborated 
by the real-space observation of SkX in Feo.5Coo.5Si thin 
film in a wide magnetic field and temperature range @. 
Lorentz force microscopy showed that Skyrmions form a 
triangular lattice with the magnetization vector antipar- 
allel to the applied magnetic field in the Skyrmion center 
and parallel at the periphery, as was also concluded from 
the neutron experiment @. 

The next important step is to explore dynamics of 
Skyrmion crystals and the ways to control them in anal- 
ogy to the actively studied current- and field-driven mo- 
tion of ferromagnetic domain walls [9]. Recent observa- 
tion of the rotational motion of the SkX in MnSi suggests 
that Skyrmions can be manipulated by much smaller cur- 
rents than domain walls 10]. 

In this Letter we study the coupled dynamics of spins 
and charges in the SkX, focusing on effects of the non- 



trivial Skyrmion topology and effective gauge fields in- 
duced by the adiabatic motion of electrons in the SkX. 
We derive equation of motion for the collective vari- 
ables describing the SkX, calculate its phonon disper- 
sion, and discuss a new form of damping, which can be 
the dominant spin-relaxation mechanism in half-metals. 
In addition, we consider new transport phenomena, such 
as the topological Hall effect in a sliding SkX and the 
Skyrmion Hall effect. We also discuss the Skyrmion pin- 
ning by charged impurities and estimate the critical cur- 
rent above which the SkX begins to slide. 

Low- energy excitations in Skyrmion crystal: An iso- 
lated Skyrmion has two zero modes corresponding to 
translations along the x and y directions. Since an ap- 
plied magnetic field opens a gap in the continuum of spin- 
wave excitations, the low-energy magnetic modes in SkX 
are expected to be superpositions of the Skyrmion dis- 
placements, or the phonons. The phonon modes, as well 
as the coupling of Skyrmion displacements to the external 
current, can be consistently described in the framework 
of elasticity theory. 

We begin with the spin Hamiltonian Hg = 
Jd 3 x [^(Vn) 2 + £n • [V x n] - • n] , where J is 
the exchange constant and D is the DM coupling that 
stabilizes the SkX configuration n(x) in some interval of 
the magnetic field H = [ll|. We calculate the 'har- 

monic lattice energy' by considering a deformation of the 
SkX, n(x, t) — n(x — u(x, £)), where the collective coor- 
dinate u(x, t) varies slowly at the scale of the SkX lattice 
constant. The result is: 



H s = dnJ 



d 2 x 



[(V^) 2 + (V My ) 2 ], 
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a 4 is the 



where d is the film thickness and £ 
characteristic length scale of SkXJjJ, with a being 



2 



the lattice spacing. The dimensionless quantity r\ = 
J nc d 2 x (diii ■ diti) encodes the information about D 
and H, and is called shape factor in what follows. 

When an electron current is flowing through the metal- 
lic film, the conduction electrons interact with local mag- 
netic moments through the Hund's rule coupling Hr = 
—JnSip'(T ■ nip, where ip is the electron operator. In 
the case of small current density and the Skyrmion size 
much larger than the Fermi wavelength of conduction 
electrons, one can apply the adiabatic approximation in 
which the electron spins align perfectly with the local 
moment, ip is projected into the fully polarized state by 
ip = x|n) with a ■ n|n) =|n). Then the electron action 
£?ei = / dtd 3 x{ihip' i ip + ^ip^V 2 ip + J n Sip*(T ■ nip] can be 
rewritten as S e i — J dtd 3 x\iUx^X — eao — 27^X^(~ — 
fa) 2 x + JnSx^x], where a /1 =|f(l - cos 0)0^ with 9 
and ip being the spherical angles describing the direc- 
tion of the local magnetization [H, lla ]. The gauge po- 
tential gives rise to internal electric and magnetic 
fields, e and h, acting on spin-polarized electrons pass- 
ing through the SkX in analogy with the electromag- 
netic gauge field. Crucially, the internal magnetic field 
b = Vxa = I s (n • d x n x d y n) z is intimately re- 
lated to the topological charge Q of Skyrmions byfUJ 
Q = j- J d 2 x (n • d x n x d y n) = ±1, where the integra- 
tion goes over the unit cell of the SkX. In the language 
of internal gauge field, this topological feature is nothing 
but the quantization of internal flux Q — J h ■ dS in units 
of hc/e. The coupling of the electric current to the inter- 
nal gauge field induced by the SkX, H lnt = — - / d 3 xj ■ a, 
has a simple form in terms of the collective coordinates 
introduced above: 



HQ [<Px 

-Hint = d— J -^-(U X Jy - UyJ X 



(2) 



The crucial difference between the SkX and a con- 
ventional crystal is the form of kinetic energy. The 
spin dynamics originates from the Berry phase action, 
<Sbp = ~ Jdtd 2 x(cos9 - l)<p. Here, 7 = h{s a +x/2) 1 where 
x is the filling of the conduction band, and S + x/2 is the 
total spin averagely per lattice site. In terms of u the 
kinetic energy has the form 



Two consequences follow immediately. First, the dis- 
persion of phonons in the SkX obtained from Eq. (J4j) is 
quadratic, 



Hlo 



rjja 2 



(5) 



in contrast to the linear phonon dispersion in usual crys- 
tals and similar to the dispersion of magnons in a uniform 
ferromagnet. Since u x and u y play the role of the coor- 
dinate and momentum, the longitudinal and transverse 
phonon modes in the SkX merge into a single mode cor- 
responding to the rotational motion of Skyrmions, which 
leads to the quadratic dispersion. Secondly, the SkX can 
move as a whole driven by the charge current j, with a 
velocity V|| = u = — ^pj. This rigid motion of SkX leads 
to several interesting results discussed below. 

Hall effect due to SkX motion: In such nontrivial spin 
textures, the external magnetic field (less than 0.2T for 
MnSi) is more than one order of magnitude smaller than 
the internal one, so that it would be neglected in what 
follows. As can be seen from Eq. the collective 

coordinates u x and u y play the role of electromagnetic 
gauge potentials A y and — A x , respectively. It is thus 
expected that the temporal variation of u induced by 
the current leads to a transverse potential drop. This 
Hall-type effect can also be intuitively understood using 
the internal magnetic field b introduced above. A moving 
spin texture n(x — V||t) induces an internal electric field 
e analogous to the electric field of a moving magnetic 
flux and related to the internal magnetic field by e = 
— - [Vii x b] . For SkX with b = b z z, this electric field 
generates an electric current in the direction transverse 
to Vii resulting in the Hall conductivity: 



Act™ 



e(b z 



2S + x mc 
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Sep = — dt—(u x u y - u y u x ). (3) by b via e = \ [v x b], and a^° p /a, 

T J £ where v is the electron velocitv. Our 



This form of the Berry phase shows that the collec- 
tive variables u x and u y describing local displacements 
of Skyrmions form a pair of canonical conjugate vari- 
ables, replacing cos (9 and tp. This characteristic prop- 
erty of SkX leads to several unusual responses to ap- 
plied electric currents and fields. It originates from 
the Skyrmion topology and distinguishes SkX from non- 
topological spin textures such as spirals and domain wall 
arrays. 

Using Eqs.([T|), and ([3]), we obtain equation of mo- 
tion for u: 



where m is the electron mass and r is the relaxation 
time. The average internal magnetic field is (b z ) — 
where $o is the elementary flux and 27r£ 2 is the area of 
the unit cell of the SkX. This Hall conductivity has the 
same order of magnitude as the one resulting from the so- 
called topological Hall effect observed in a static SkX[15|. 
The latter effect is nothing but the Hall effect induced 

s e(6 2 )r/mc, 

where v is the electron velocity. Our new effect differs 
by the factor of — 2 g +x from the topological Hall effect. 
Its physical origin can be easily understood by noting 
the total force acting on a single conduction electron is 
F = — -[(v — Vy) x b], i.e. the Lorentz force on elec- 
trons due to the internal magnetic field of the SkX de- 
pends on the relative velocity of electrons and Skyrmions. 
When the SkX begins to slide above the threshold elec- 
tric current j c [lB|; the net topological Hall voltage will 
be suddenly reduced by the factor , which is how 
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Q^i x V 2 u. 
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(4) 



2S+a;' 

the effect of the spin-motive force and the collective shift 
of Skyrmions can be identified experimentally. 

New damping mechanism and Skyrmion Hall effect: 
Previously we have systematically discussed the novel 
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effects related to the internal magnetic field. A natu- 
ral question thus arises as to whether there is any new 
phenomena associated with the intrinsic internal electric 
field, which is e$ = — diao — -Cbi = £ (n • c^n x ri). Due 
to the time derivative in this expression, its effect is ab- 
sent in the static spin texture. However, in the present 
case, the motion of SkX makes it nonvanishing, and leads 
to an additional current j' by j' = ere with a the con- 
ductivity of electrons. Substituting this current into the 
Landau-Lifshitz-Gilbert equation [l2L [l3| 



SH S 



Sn 



• a n x n 



(7) 



the time derivative n receives a correction given by 

Sh = (e • V)n = a' (n • diin x n) c^n. (8) 

The corresponding dimensionless damping constant is 
a> = (2S+x) etfag^c ' wnere a fs ~ 1/137 is the fine struc- 
ture constant. The time derivative in the r.h.s. of Eq.© 
shows that the current induced by internal electric field 
leads to dissipation. In contrast to Gilbert damping this 
new mechanism does not require relativistic effects and 
only involves the Hund's rule coupling that conserves the 
total spin. The relaxation of the uniform magnetization, 
described by Gilbert damping, is clearly impossible with- 
out the spin-orbit coupling, which breaks the conserva- 
tion of the total spin|17jj. This argument, however, does 
not apply to inhomogeneous magnetic textures where the 
breaking of the rotational symmetry by noncollinear spin 
orders enables the relaxation without the spin-orbit cou- 
pling (note that a' vanishes as £ — > oo). Despite the 
non-relativistic origin, a' depends on the DM coupling, 
as the latter determines the Skyrmion size. Estimates of 
a' made below show that in half-metals it can greatly 
exceed a. 

The effect of this new dissipation can be observed by 
tracing the trajectory of Skyrmion motion. Including the 
new dissipation term, the modified equation of motion Q 
for the rigid coilective coordinates u(t) has the form 

u = -— j - Q (ar) + a'rf)i x u, (9) 

where the second shape factor vj is given by rj = 
& J uc d 2 x (n • d x n x d y n) (d.n ■ <9;n) / J uc d 2 x (<%n ■ 9;n). 
The new dissipation term in Eq.([9]) is obtained by mul- 
tiplying Eq.© with djn, using ri = — (u • V)n, and 
integrating over one unit cell. The whole damping term 
leads to a transverse motion with velocity 



V 



Q(aT] + a'r)') 



[V|| 



(10) 



This Skyrmion Hall effect can be observed by real-space 
images of Lorentz force microscopy. The corresponding 
Hall angle is 9 — arctan(ary + ot'rf). The estimate given 
below shows that main contribution to 6 comes from the 
new dissipation mechanism. 



Pinning of Skyrmion crystal: Next we consider the 
pinning of the SkX by charged impurities. The pinning 
results from spatial fluctuations of the impurity density 
and variations of the spin direction in the SkX. Variations 
of the density of charged impurities Snt give rise to local 
variations of the electron density n e and since the double 
exchange constant J is proportional to the latter, we have 
5J ~ J5ni/n e . The energy per Skyrmion Eg ~ Jd/a. 
Denote the number of impurities in this volume by N\ 
with (Ni) — Uil-n^d and the variance 6Ni = y 



/iVi, we 

obtain the typical variation of the Skyrmion energy: 



Vi 



SJ d - 

a 



J 




n e 2n^ 2 a n e at; 

The potential energy density is then V = Vi/(27r£ 2 ). 
Substituting n,; ~ (la ) , where I is the electron mean 

free path, and n e = Jr, we obtain V ~ ^yi^x \ft 
The pinning regime of the whole SkX depends on the ra- 
tio of the pinning energy V\ and the elastic energy Eg of a 
single Skyrmion 18]. Let L 2 be the number of Skyrmions 
in the domain where u ~ ^. The energy gain due to 
the impurity pinning in the domain is ~ —V\L, while 
the elastic energy cost ~ — is independent of the do- 
main size. Minimizing the total energy per Skyrmion, 
~ i, we obtain L ~ L>1 corresponds to the 
case of weak (or collective) pinning of SkX, while L ~ 1 
corresponds to the strong pinning regime. 

The pinning potential gives rise to the spin transfer 

torque — Q%d x f^] m tne right-hand side of Eq.(gJ|. 
In the steady state of moving SkX this torque has to be 
compensated by the interaction with the electric current. 
The critical current density is then 



h d 



dV 
du 



h dL ' 



3c~T-{ — ) -T — ■ (I 2 ) 

steady state 

in the weak pinning regime, while in the strong pinning 
case L has to be substituted by 1. Similarly, one can 
estimate the gap in the spin wave spectrum due to the 
pinning: 



h^ 2 /d 2 V\ h^ 2 V L a 3 V 



pin 



du 2 



d L 2 £ 2 dS L 



(13) 



Estimates: For estimates we consider MnSi where Mn 
ions form a (distorted) cubic sublattice with a = 2.9 A. 
The length of the reciprocal lattice vectors of the SkX 
- 0.035 A" 1 corresponds to £ ~ 77 A and D - 0.1J. The 
kinetic energy scales as Ti 2 /m£ 2 < Jh ~ leV, so that the 
adiabatic approximation is justified. The electron density 
n e — 3.8 ■ 10 22 cm _3 [l5j corresponds to x = n e a 3 ~ 0.9 
charge carriers per lattice site, while the residual resis- 
tivity p ~ 2[iQ ■ cm[l5j gives an estimate of the impurity 
concentration Xi ~ 5- 10 -3 . From the magnon dispersion 
in the spiral state[ll|, J ~ 3 meV. 

Using these parameters we get a' ~ 0.1, which shows 
that the damping resulting from the electric currents gen- 
erated by non-collinear spin textures can be the domi- 
nant mechanism of spin relaxation in half-metals, where 
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the Gilbert damping constant a is one-three orders of 
magnitude smaller [20. l2lj|. We note, however, that since 
the typical electron mean free path £ ~ 500 A is larger 
than the Skyrmion size £, the relation between the cur- 
rent and internal electric field is nonlocal. The topo- 
logical Hall angle, H ~ = — L -?r £ at T = 
and the change in #h induced by the sliding Skyrmion 
crystal [see Eq.©] are also ~ 0.1. For a lOnm thick 

film the parameter L ~ x ^^ 2 ? dl ~ 10 3 , i.e. the pin- 
ning of Skyrmions by charged impurities is exceedingly 
weak and the corresponding values of the pinning fre- 
quency huip\ n ~ 5 • 10~ u meV and the critical current 
j c ~ 0.2A • cm -2 . This value is much lower than that 
for domain walls Q and smaller than the ultralow thresh- 
old current j c ~ 10 2 A • cm -2 observed in bulk MnSiflol]. 
which may be attributed to other pinning mechanisms 
and the different dimensionality of the system. This low 
current density also justifies the adiabatic approximation 
applied in this work. 

Comparison with vortex dynamics: Finally, we com- 
pare the dynamics of SkX with that of the vortex lines 
(VL) in type II superconductors. A similar quadratic dis- 
persion was obtained for VL 22]. However, in supercon- 
ductors it results from long-ranged interactions between 
vortices, while the interactions between Skyrmions are 
short-ranged (if one ignores the relatively weak dipole- 
dipole interactions). The absence of long-range interac- 
tions in SkX ensures the stability of the quadratic disper- 
sion. Furthermore, the kinetic terms in VL and SkX are 
completely different. For VL u x and u y are two indepen- 
dent variables, while for SkX the are conjugated variables 



as in Eq.Q. Therefore VL are massive, while Skymions 
are not. When a supercurrent flows through the type II 
superconductor, the charged Cooper pairs are deflected 
by the VL through the Lorentz force, which in turn gives 
rise to the transverse motion of the VL. The VL dynamics 
is usually assumed to be overdamped[23j, the kinetic en- 
ergy of VL is neglected, and the Lorentz force is assumed 
to be counterbalanced by the friction force. In contrast, 
the damping of Skyrmions is relatively weak. The spin 
torque resulting from the strong Hund's rule coupling re- 
sults in a nearly longitudinal motion Skyrmion motion. 
The Hall motion of VL results in a longitudinal voltage 
drop, which is not important for SkX motion due to the 
small Hall angle. 

During the completion of this paper we became aware 
of a recent paper by Kim and Onoda addressed the dy- 
namics of Skyrmions in an itinerant double-exchange fer- 
romagnet using a Chern-Simons-Maxwell approach (2^. 
Their focus, however, seems to differ from ours. We are 
grateful for the insightful discussions with Prof. Yoshi- 
nori Tokura. This work is supported by Grant-in- Aids for 
Scientific Research (No. 17105002, 19019004, 19048008, 
19048015, and 21244053) from the Ministry of Education, 
Culture, Sports, Science and Technology of Japan, and 
also by Funding Program for World-Leading Innovative 
R&D on Science and Technology (FIRST Program). JZ 
is supported by Fudan Research Program on Postgrad- 
uates. MM is supported by the Stichting voor Funda- 
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